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I. INTRODUCTION 

The study of the optical properties of nanoscale semi- 
conductor structures, such as small clusters, nanocrys- 
tals or polymer chains, is acquiring increasing importance 
also in view of optoelectronic applications. In these sys- 
tems, however, the theoretical description of the opti- 
cal response faces additional difficulties with respect to 
extended crystals, because correlation effects are known 
to be enhanced. Besides the single-particle energies and 
wave-functions for electrons and holes, also an appro- 
priate description of Coulomb interactions between elec- 
trons and holes is required. For this reason semiempiri- 
cal or approximate methods (e.g., based on the effective- 
mass approximation) were often the only feasible treat- 
ment so far. In many cases of interest, however, such 
methods are not sufficiently accurate and predictive. 

A strong impulse towards ab-initio approaches is re- 
cently coming from the success of standard Density Func- 
tional Theory (DFTV,methods to treat extended three- 
dimensional systems,™ with unit cells that may now in- 
clude thousands of atoms. In semiconductors, however, 
such DFT approaches are known to account only insuffi- 
ciently for Coulomb correlation effects. The main reason 
for this shortcoming is that the dynamics of electron- 
hole excitations is not only governed by the DFT single- 
particle states but also by polarization effects, which are 
not fully accounted for in a standard DFT calculation. 
More specifically, in semiconductor carriers polarize their 
surrounding; in turn, this induced polarization affects 
both the propagation of carriers (quasiparticlc rcnormal- 
izations) and the Coulomb attraction between optically 
excited electrons and holes (dielectric screening). 

Within the framework of Green's functions it has been 
known for a long time how to correct for such polar- 
ization screening. For instance, Hedin and coworker (see 
Refs. |[||for a review) showed how to consistently develop 
a set of equations accounting for Coulomb correlation ef- 
fects. In lowest order, this set of equations reduces to 
the calculation of the screened Coulomb potential W, ob- 
tained from the solution of Dyson's equation for W within 



the random-phase approximation, and the self-energy S 
which is calculated within the so-called GW approxi- 
mation (GWA).Btl While the GWA approach provides 
a strongly improved description of single electrons and 
holes, the description of optical absorption requires in 
addition a proper treatment of the attractive Coulomb in- 
teraction between the optically excited electron-hole pair. 
This electron-hole correlation is conveniently accounted 
for by solving the Bethe-Salpeter equation, which de- 
scribes the propagation the two-particle Green's function 
(for details see Ref. || and references therein). 

On the other hand, within the last couple of years 
phenomenological semiconductor optics, i.e., the simpli- 
fied bandstructure. and dielectric screening-description of 
semiconductors, B'B has become an extremely vivid and 
successful field of research. Here, the main interest is on 
the description of the nonlinear optical response of semi- 
conductors and semiconductor nanostructures, and of op- 
tically induced coherence effects. A particularly trans- 
parent description scheme for such systems is provided by 
the dynamics-controlled truncation of the density-matrix 
hierarchy^ Contrary to the Green's function approach, 
which describes by means of many-body perturbation 
theory the dynamic propagation of quasiparticles, within 
the framework of density matrices one selects an appro- 
priately chosen set of one-time correlation functions. If 
Coulomb interactions conserve the number of electron- 
hole pairs, it becomes possible to rigorously classify these 
correlation functions in powers of the electric field of the 
exciting laser,B thus providing a consistent description 
scheme for the optical properties of semiconductors. 

If, instead of adopting from the beginning a simplified 
semiconductor description, we start from the results of a 
first-principles calculation, the most interesting question 
arises: How is it possible to describe dielectric screening 
and Coulomb correlation effects within a first-principles 
density-matrix description? To the best of our knowl- 
edge, up to now this problem has not been addressed in 
the literature. Its answer, however, would be very inter- 
esting for a number of reasons. First, it would be helpful 
to establish a link between the fields of semiconductor op- 
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tics and of first-principles calculations; more specifically, 
such synthesis would allow the first-principles calcula- 
tion of non-linear optical properties of semiconductors. 
Second, the framework of density matrices requires only 
some basic knowledge of quantum mechanics and statis- 
tical physics — in contrast to the framework of Green's 
functions, which rests on a large machinery of defini- 
tions and auxiliary functions. Thus, one might hope that 
certain types of approximations become more transpar- 
ent; furthermore, the respective strengths of the two ap- 
proaches might be helpful in establishing approximation 
schemes beyond the presently used approaches. 



This paper is devoted to the discussion of a first- 
principles calculation of dielectric screening within the 
framework of density matrices. It will be shown that 
non-pair conserving Coulomb couplings, whose neglect is 
at the heart of the dynamics-controlled truncation ap- 
proach, have to be retained for a proper first-principles 
density-matrix description of dielectric screening and 
quasiparticle renormalizations in optically excited semi- 
conductors. Adopting the usual random-phase approxi- 
mation and keeping our analysis at the lowest level of ap- 
proximation, we will propose a method for the ab-initio 
calculation of optical absorption which closely resembles 
that of a corresponding Green's function analysis: first, 
the exciton energies E x and wavefunctions ^S> x are ob- 
tained from the solutions of the two-particle electron-hole 
Schrodinger equation, where the electron and hole inter- 
act with each other through the screened direct Coulomb 
term and the unscreened exchange one; second, the op- 
tical transition energies are given by E x whose oscillator 
strengths are calculated from ^> x . As the main difference 
with respect to the corresponding Green's function anal- 
ysis, our prescription suggests to use static screening in 
the calculation of the direct Coulomb term. 



We have organized this paper as follows. In Sec. II we 
present our model Hamiltonian and introduce the frame- 
work of density matrices; a short overview over one of the 
most widely used truncation schemes in phenomenologi- 
cal semiconductor optics, the dynamics-controlled trun- 
cation, will be given. Sec. Ill is devoted to the analysis 
of dielectric screening within the framework of density 
matrices; the pertinent dynamic variables and approxi- 
mation schemes are identified. To keep our analysis as 
simple as possible, details of our calculations are post- 
poned to Appendices A and B. Finally, in Sec. IV we 
compare our results with those obtained within the GWA 
and Bcthc-Salpeter-equation approach and draw some 
conclusions. Finally, in Sec. V we draw some conclu- 
sions and summarize our prescription for the ab-initio 
calculation of optical absorption in semiconductors. 



II. DENSITY MATRICES 
A. Hamiltonian 

Let us assume that the single-particle energies e a and 
wavefunctions 4> a (f) of the semiconductor system under 
investigation are obtained from a DFT calculation. Al- 
though we shall not be too specific about the details of 
the underlying ab-initio approach, we assume that the 
Hartree and exchange interactions are included without 
approximations^ a discussion of differences with respect 
to the local-density approximation (LDA) will be given 
below. Let H Q denote that part of the Hamiltonian which 
includes the kinetic energy and the lattice interactions, 
and H' that part of the electron-electron Coulomb in- 
teractions H' which are included through the Hartree, 
exchange, and correlation potentials vh -\-v x + v c . Then, 
e a and (jfaif) are the DFT eigenstates of H a — H Q + H', 
and SH = H' — H' are those Coulomb interactions which 
are not included in the DFT approach (and which will be 
treated in this paper by means of perturbation theory) . 
The basic Hamiltonian describing the carrier system is 
thus of the form: 



H = H + H op + SH 



(1) 



with the single-particle Hamiltonian H Q \ the light- 
semiconductor coupling -Hop! and the Coulomb interac- 
tions SH which are not included in DFT. Note that for 
simplicity in Eq. (|l|) we have neglected all types of carrier- 
phonon interactions or other environment couplings. In 
the spirit of the GWA approach, in this work we consider 
in the subtractions to SH only the Hartree and exchange 
terms, thus ignoring possible double counting of Coulomb 
correlations v c . Hence, with the Fermionic field operators 
i/>(r), which create an electron at position r, we obtain: 



5H ~ \ J d ( rr 



, ^{r)ip\r')xp(r')ip(r) 



d{rr') [vjj(r)6(r — r 1 ) + v x (r, r')]xp* (r')ip(r). 

(2) 



In the DFT semiconductor groundstate all states be- 
low the fundamental bandgap are occupied (valence-band 
states) and all states above the bandgap are unoccupied 
(conduction-band states). Next, we transform to the 
electron-hole picture and introduce the Fermionic field 
operators and , where: c\ creates an electron in the 
conduction-band state 4>\{r) with energy e\ (henceforth 
we use numbers 1,2,... to label single-particle states); 
d\ creates a hole in the valence- band state (fi^{r) with 
energy e\. Since the hole describes the properties of 
a missing electron in the valence bands, we adopt the 
usual definitions that for a valence band state v the cor- 
responding hole state 2 is related through e|} = — £ v and 
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( t ) 2 ( r ) = [ < / ) i'(' r )]* (with a properly chosen zero-point en- 
ergy). With these operators: the single-particle Hamil- 

tonian is of the form H a — J2i e i c \ c i + S2 e 2 d\d 2 ; 
the light-matter couplinfi-within the dipole and rotating- 
wave approximations is:Ell3 

H o P = -\£o M 2i d 2 c, + e-*-* 44), (3) 

12 

with £ exp ±ia;i the electric field of the exciting laser 
light and Mfij the optical dipole matrix elements; the 
Coulomb terms can be split into three parts SH = 
5HW + SH^ + SH^fi where SH^ conserves the num- 
ber of electron-hole pairs, SH^ 1 ' changes the number of 
electron-hole pairs by one, and 5H^ 2 ' by two (see Ap- 
pendix |A| for details). 



and screens the Coulomb interactions in SH^ by the 
static dielectric constant (see Sec. Ell for a first-principles 
motivation of such approximation). Hence, the remain- 
ing terms H a + SH^ conserve the number of electrons 
and holes, and the only source for the creation or destruc- 
tion of electron-hole pairs is through the light-matter cou- 
pling H op ; more specifically, inspection of Eq. (^|) reveals 
that the creation (destruction) of electron-hole pairs is 
through terms of the form £ or t^d' i (£ a dc). In the pio- 
neering work of Axt and StahO (see Ref. || for a review) 
the authors first noted that, since any pair of field opera- 
tors c' d) or dc comes with an electric field £ , the density 
matrices can be classified according to their power in the 
electric field. 

Thus, in linear optical response there is only one cor- 
relation function 



B. Density- matrix hierarchy 



Y 12 = (d 2Cl ) 



(5) 



Having established our model in Sec. II A, we next 
discuss how to treat Coulomb correlation effects due 
to SH. Within the framework of density matrices the 
central quantities are the one-time correlation functions 
(A)t — tr(p t _4), with p t the statistical operator at time 
t, tr() denoting the trace over a complete set of ba- 
sis functions, and A an operator consisting of c, d, c\ 
and d); for instance, {c s c) t {{d^d) t ) describes the occu- 
pation of electron (hole) single-particle states at time 
t, whereas higher-order correlation functions, such as, 
e.g., (c ! a ! dc)t, account for correlations among two or 
more carriers. As the central approximation within any 
density-matrix description one has to restrict oneself to 
a limited number of correlation functions. The temporal 
evolution of the correlation ftmctions {A) t is then pro- 
vided by Ehrenfest's theorem! 



d t (A) t = ([A,H]) t) 



(4) 



which, together with the restriction to a finite number 
of correlation functions, provides the basis of the frame- 
work of density matrices. Before addressing the prob- 
lem of dielectric screening in Sec. HI, in the following we 
briefly review one of the most commonly used trunca- 
tion schemes in phenomenological semiconductor optics: 
The dynamics- controlled truncations^ This discussion will 
help us to establish some of the basic concepts and nota- 
tions. 



C. Dynamics-controlled truncation 

In describing the optical properties of conventional 
semiconductors one often relies on effective models. Be- 
sides a simplified bandstructure description in terms of 
fe.p-tjieory or effective-mass approximations one usu- 
ally :tlu employs the envelope- function approximation; ne- 
glects the non-pair conserving terms SH^ and SH 1 * 2 " 1 ; 



which is of lowest order in £ a . Y describes the correlation 
between an optically excited electron and hole (note that 
Y is a non-equilibrium quantity which is nonzero only 
in presence of an exciting laser and which vanishes in 
thermal equilibrium). Its relation to the total interband 
polarization V{t) is given through V(t) = ^2 12 M|ili2,i, 
with M oh the optical dipole matrix elements. Recalling 
that the optical absorption is related through Maxwell's 
equation to the total interband polarization ^(i),™ one 
immediately notes that the microscopic interband polar- 
izations Y completely determine the linear optical re- 
sponse of semiconductors. Their temporal evolution is 
given bjt,Eq. (jl]), and we find after some straightforward 
algebra:ElEI 



ilia - (e? + e\)Y 12 - \£ Q e 



eh 



fV ee,hh - V° h ' he ) Yin 
V v 11,22 v 12,21'' 12 ' 



(6) 



where we have implicitly assumed summation over barred 
indices 1 and 2. On the right-hand side of Eq. (^): the 
first term accounts for the single-particle states avail- 
able for adding an electron-hole pair; the second term 
describes the creation of electron-hole pairs through cou- 
pling to the light field; finally, the third term accounts 
for the electron-hole Coulomb interactions (with V ee ' hh 
the direct and V eh ' he the exchange term), and is respon- 
sible for excitonic renormalizations in the optical spec- 
tra (which would be absent in a simple-minded single- 
particle description). 

Eq. (H) is conveniently solved by finding the polariza- 
tion eigenmodes. The homogeneous part of Eq. (|) is 
then transformed to .an-eigenvalue problem ("excitonic 
eigenvalue problem" ) ,ErEla whose solutions provide the ex- 
citon energies E x and wavefunctions ^f 2 - In the optical 
spectra the optical transition energies are then provided 
by E x , and the respective oscillator strengths are given 
by IE12 *12 M 2i! 2 - 
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III. DIELECTRIC SCREENING 

As it is well known that dielectric screening is of cen- 
tral importance for the quantitative description of the 
optical properties of semiconductors, within an ab-initio 
approach it would make little sense to directly employ 
Eq. (^) with the bare (i.e., unscreened) Coulomb matrix 
elements V. For instance, in conventional semiconduc- 
tors, such as Si or GaAs, dielectric screening leads to a 
reduction of the unscreened V's by a factor of approxi- 
mately ten. We are thus faced with the central question 
of this work: How can dielectric screening be described 
within the framework of density matrices? 

Quite generally, dielectric screening is a process where 
a carrier polarizes its surrounding medium. In turn, a 
second carrier not only feels the bare Coulomb interac- 
tion exerted by the first carrier, but also the field pro- 
duced by this induced polarization cloud. In an electron- 
hole picture, such polarization effects are described by 
(virtual) excitations of electron-hole pairs which result 
in microscopic polarization fields. Thus, a proper first- 
principles description of dielectric screening requires be- 
sides the pair-conserving Coulomb terms 6H(°> also the 
non-pair conserving terms <5ijW and 5H {2 \ 

Inclusion of such terms, however, spoils the concept of 
the dynamics-controlled truncation of the density-matrix 
hierarchy, since electron-hole pairs are no longer solely 
created (destroyed) by the light field but, in addition, 
also by Coulomb interactions. Thus, it is no longer pos- 
sible to rigorously classify correlation functions according 
to their power in the electric field £ a . However, since both 
the light field and the Coulomb interactions create elec- 
trons and holes only pairwise, we can still assume that 
the classification of correlation functions in numbers of 
electron- hole pairs provides a meaningful concept. For 
simplicity, in this paper we restrict ourselves to the case 
of correlation functions with at most two electron-hole 
pairs, and we show that already at this level of approxi- 
mation dielectric screening is included. 

A. Dynamic variables 

Besides the interband polarizations Y12, we hence need 
the electron (hole) distribution function C12 = (<4ci) 

{D12 rp (^2^1)) an d the two-particle correlation func- 
tions 

■Bl2,34 = ( d 4 C 3 d 2 C l) ( 7 ) 
^12,34 = (44<W. (8) 

In analogy to the physical meaning of Y as a measure of 
the electron-hole correlation and C (D) as the occupation 
and intraband coherence of single-particle electron (hole) 
states, we can approximately interpret B as a measure 
for the coherence between two electron-hole pairs and N 
as the occupation of electron-hole pairs. 



Next, we derive the equations of motion for the various 
correlation functions. Before doing so, we employ a fur- 
ther approximation. From the analysis of the nonlinear 
coherent optical response (i.e., the dynamics-controlled 
truncation at the level of which involves Y, C, D, B, 
and N)U it is known that the main effect of the carrier 
distribution functions C and D is to reduce the optical 
transition rates because of state filling (i.e., Pauli block- 
ing) . As we expect such effects to be of minor importance 
for the problem of our present concern, in the following 
we shall neglect C and D. Thus, our set of dynamic 
variables is provided by Y, B, and N. 

B. Equations of motion for Y, B, and iV 

The program pursued in the following is the deriva- 
tion of the equations of motion for Y, B, and N. This 
is done by using Ehrenfest's theorem, Eq. (|J), together 
with the definitions for Y, B, and N, Eqs. (^,0,||) (see 
also Appepdix |a|). After some straightforward algebra 
we obtainO 

iY 12 S (e- + 4) Y 12 - \E e- wt M c £ 

_/yCO,hh _ ych.hcx Y 
V v 11,22 v 12,21/ 12 

. /Tree, eh iree.ch^ , ycc,hc R 

+ ( - v ii,34 - V 31,14^ JV 12,34 + V 11,43 ^12,34 

_ ( V 22,34 ~ V 42,32) ^12,34 - V 22\43 B 12,34 (9) 

(here and henceforth we implicitly assume summation 
over all single-particle indices with exception of those ap- 
pearing on the left-hand side) . Comparison with Eq. (|g) , 
which was derived by neglecting non-pair conserving 
Coulomb couplings, reveals the appearance of an addi- 
tional coupling to the two-particle correlations B and 
N, which is mediated through the non-pair-conserving 
Coulomb interactions SH^ and SH^ 2 \ As will be shown 
in the following, such terms are needed for a proper de- 
scription of dielectric screening. 

Let us discuss this coupling between Y and B, N in 
slightly more detail. Suppose that the system under in- 
vestigation is a bulk semiconductor. Hence, the single- 
particle states 1,2,... consist of a wavevetor k and a 
band index n. As the light field couples electron and 
hole states with approximately opposite fc-values, the 
optically induced interband polarization is of the form 
Yk = {d_ k c k }, where for notational simplicity we have 
dropped all band indices. The coupling VB in Eq. (|^) 

is then of the form E P q v (<?)( d - P +g c P d -k+ q c k)^ with 
a corresponding expression for VN. Hence, these terms 
describe how the propagation of Y is modified by the 
presence of polarization fluctuations in the system; as we 
will show in the following, such fluctuations are induced 
by the interband polarization Y, thus resulting in a self- 
interaction-like process where: the interband polarization 
Y induces polarization fluctuations through the non-pair- 
conserving Coulomb couplings SH^- 1 '; these fluctuations 
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propagate in time, which is described by the equations 
of motion for B and N (to be derived below) ; and finally 
couple back to Y, which is described by the terms VB 
and VN in Eq. (|). 

Next, we derive the equations of motion for B and N. 
To keep our analysis as simple as possible, we employ 
from the beginning a number of well-controlled approxi- 
mations. First, we assume that before arrival of the excit- 
ing laser no two-particle correlations B and N are present 
(as will be discussed below this is only an approximation) . 
Hence, Y, B, N are induced by the light-semiconductor 
coupling H op , and in linear response it suffices to keep in 
the dynamics only terms linear in Y, B, N (i.e., we ne- 
glect contributions proportional to, e.g., Y 2 or YB). Sec- 
ond, we neglect the light coupling in the equations of mo- 
tion for B and N (such terms would describe interference 



terms between dielectric screening and light coupling). 
Third, we shall emj>lojy-|the so-called random-phase ap- 
proximation (RPA):tL3li3 For .81,2,34 an d iVi.2,34 we keep 
only Coulomb terms which scatter one pair of particles 
12 (34) and leave the other pair 34 (12) unaffected. In 
other words, we assume that polarization fluctuations 
created through light-coupling and Coulomb processes, 
respectively, move independently of each other. In a 
fc-space representation for the polarization fluctuations 

( d -P+g c P d -k+ q c k) and ( c i d -p+ g d -k+ q c k) one readily 
observes that these RPA terms contain all those contribu- 
tions which involve Coulomb matrix elements V(q) with 
a momentum exchange of q which is independent of k 
and 

With these approximations we then find: 



iB l2M S (e\ + e\ + e§ + e\) B l2M 

i /-i j-ch.hc -v rcc,hh\ ri I /-ch,ch -i 7-ch,ch\ * j 

+ ^ V 12,21 ~ V 11,22V ^12,34 + l V 12,12 ~ V 12,l2^ JV 34,12 

i / x rch.hc -i rcc.hhx j-y . /-v 7-ch,ch -v 7-ch.ch\ j\ T 

+ l V 34,«3 " V 33,4i) ^12,34 + 1* 3 4,U ~ V 3i,U> ^ 12 ^ 

i /-ircc,ch Trcc,ch\ ^ /Trhh,eh xrhh,eh\ \r _ 

+ ^ V ll,34 _ W 31,14> J 12 - <> V 22,34 ~~ V 42,32- ) J 12 

j Arcc,eh -»rcc,ch\ ~\/*_ /T7"hh,ch -irhh.chx y i /-v/-ch,ch -treh,eh\ f-\c\\ 

+ ^ V 33,12 ~ V 13,32^ ^34 - (V^ ^ - V 2i 14 J r 3 4 + ^^ 34 - V 14i32 j (1UJ 

iN l2M S (t\ + 4~el~e\) N 12M 

I /-i rCh.llC t rCC.hh\ \T _ i f\ rCh,ch t rCh.cllN 73* 

+ l V 12,21 ~ V 11,22V iV 12,34 + I V 12,12 _ V 12,12^ ^12,34 

- - Vg£) ^12,34 - CV^S - Vj££) 5l2,34- (ii) 



Let us analyze the various contributions to these equa- 
tions in slightly more detail. On the right-hand side of 
Eq. (pX|): the first term corresponds to the free prop- 
agation of B; the terms in the second and third line, 
respectively, describe Coulomb renormalization processes 
in the propagation of B (VB) and Coulomb couplings be- 
tween N and B (VN); the terms in the last two lines are 
the source terms which describe how polarization fluctua- 
tions are created through coupling to the interband polar- 
izations Y. In addition, there is a term (V°2 34 — Vf 4 '| 2 ) 
which describes the buildup of polarization fluctuations 
even in absence of light couplings. Quite generally, its 
appearance is not unexpected since up to now we have 
assumed that the semiconductor groundstate is obtained 
by simply filling DFT single-particle states. Acting with 
8H^ on the DFT vacuum, however, we immediately 
observe that Coulomb interactions create electron-hole 
pairs and thus lead to renormalizations of the DFT vac- 
uum. As a first approximation, in this paper we assume 
that such renormalizations are not of crucial importance 
for the description of dynamic processes in the propaga- 
tion of Y, and we thus neglect this term. 

In the equation of motion for N, Eq. (Ill]), we observe 



in analogy to the dynamics of B single-particle contri- 
butions (first line) and Coulomb renormalizations and 
couplings (second and third line), but no source terms 
are present. Similarly to the semiconductor Bloch equa- 
tions, Era which describe the non-equilibrium optical re- 
sponse at the level of two-point functions (i.e., C, D, 
and Y) and where the light field drives the electron-hole 
coherence Y which, in turn, acts as a source term for 
the electron (hole) occupations C (D), we here have the 
situation that Coulomb coupling first drives the electron- 
hole pair coherence B which, in turn, acts as source for 
the electron-hole pair occupation N. 

Eqs. (PHli"|) are the central equations of this work. 
Before showing how to solve this set of equations, we 
introduce two further approximations. First, we keep 
in Eqs. ( |Io| , ^Tl) only self-interaction like processes. The 
electron-hole coherence Yj 2 (Y 12 ) between the states 12 
(12) is initially created through the coupling to the light 
field. In a self-interaction process: this coherence is 
transferred through Coulomb coupling of the electron 
(hole) to 12 and to a second electron-hole pair 34 [fourth 
line in Eq. (^0|)]; the second electron-hole pair 34 prop- 
agates in presence of Coulomb renormalizations [third 
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lines in Eqs. ([[(j[Ll|)] and the initial pair 12 remains un- 
scattered; finally, this coherence between electron-hole 
pairs 12 and 34 affects the propagation of Y [third 
and fourth line in Eq. (H)]. In addition to these self- 
interaction like processes, there exist also scattering-like 
contributions: here cither the optically induced coher- 
ence between electron-hole pair 34 is transferred through 
Coulomb coupling to the pair 12 [fifth line in Eq. (10)] 
or the coherence between pair 12 is scattered to 12 [sec- 
ond lines in Eqs. ( |io| , pl]) 1 . Such terms, which will be ne- 
glected in the following, describe higher-order Coulomb 
terms (within the framework of Green's functions they 
would correspond to vertex corrections); a more detailed 
discussion of such scattering contributions has been given 
in Ref. [If 

Our second approximation concerns the neglect of 
electron-hole exchange interactions in the dynamics of 
B and N: in Eqs. ( |Io| , |lT| ) we only keep the first Coulomb 
terms in parentheses. One readily observes in a k- 
space representation that for the polarization fluctua- 



k+q C ki 



these 



tions (d_ p+q c p d_ k+q c k ) and 

terms correspond to Coulomb interactions V(q) where 
the exchanged momentum q is independent of the mo- 
menta k and p; the remaining terms involve exchanged 
momenta depending on k and p. The more general ex- 
pressions of Eqs. (|l(^,[ll]) which contain both the direct 
and exchange contributions have been given for the pur- 
pose of our later discussion about the screening of the 
electron-hole exchange interaction. 



Details of our solution scheme for Eqs. (|10|) and ( |11| ) 
are presented in Appendix |b| The key to the solution is 
the Dyson-like character of the equations of motion for 
B and N. In fact, solving Eqs. ( ^0|JTl| ) by iteration we 
observe that the repeated action of V(B + N) precisely 
reproduces Dyson's equation for the screened Coulomb 
potential W = V[l + VP + (VP) 2 + . . . ] = V + WPV; 
here P is the usual retarded polarization function within 
random-phase approximation. Finally, we arrive at (for 
details see Appendix fi|): 



J 



iY 12 = (ef + e\)Y l 



■(V 



ee,hh 
11,22 



12 



V 



ee.ee / 

n,i>i 



12 



Yi. 



V^i^Y^-^e-^M^ 



12 



U 22 ( W 12)+^22 ,11 ("15) 



-V 



eh,he\ 
12,21/ 



(12) 



with 



IV. DISCUSSION 



V(r,r';u;) 



d(ss'Q) W(r, S ; W )^^^V(s',r'), 



u — uj + i0+ 



(13) 



[the relation between the real-space and sing le-p article 
representations of V is in analogy to Eq. (A4)]; D 
is related to the imaginary part of the polarization 
function through D (r,r';uj) = — 7r -1 3P(r, r'\ u> > 0) 
[Eqs. (B5,B8)]. A particularly simple expression of 
Eq. (|l2) follows for static screening, i.e., for: 



V(r,r';cj = 0) = |[W(r, r'; 0) - V(r, r% 



(14) 



where the expression in the second line of Eq. (U3) re- 
duces to [Wg(0) - V^Jj] F I5 . Thus, in Eq. Q the 
D-terms in the first line describe renormalization effects 
in the propagation of single electrons and holes, respec- 
tively, and the terms in the second line account for the 
Coulomb coupling between the optically excited electron- 
hole pair. 



Equation (|12|) is our final expression. Before analyz- 
ing it in more detail, let us briefly recall the approxima- 
tions that were adopted in its derivation. Firat, within 
the spirit of the dynamics-controlled truncationQ only dy- 
namic variables accounting at most for two electron-hole 
pairs were considered. The electron and hole distribu- 
tion functions C and D were neglected, since from related 
work it is known that their main influence is the blocking 
of optical transitions due to state filling. Thus, our set of 
dynamic variables consists of the interband polarizations 
V, and the two-particle correlations B and N which ac- 
count for the coherence between electron-hole pairs and 
their occupation, respectively. We assumed that before 
arrival of the exciting laser light V, B, and N vanish. In 
the equations of motion for £L and N: we adopted the 
random-phase approximatior£3 which assumes that po- 
larization fluctuations propagate independently of each 
other; we neglected the direct coupling of B and N to 
the light field; only self-interaction-like processes were 
kept; finally, exchange-type interactions were neglected. 

Quite generally, none of these approximations is com- 
pulsory, and all of them could be lifted without introduc- 
ing major conceptual modifications. However, we have 
kept our analysis as simple as possible in order to em- 
phasize the two primary goals of this work: first, the 
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identification of those terms and approximations which 
are required within a density-matrix formalism in order 
to describe dielectric screening; and second, to discuss 
the respective differences between our approach and the 
combined GWA and Bethe-Salpeter equation approach. 

Within our density-matrix description, the central 
quantities are the optically induced interband polariza- 
tions Y. In general, their knowledge allows the full 
calculation of the op tical absorption spectra. However, 
as discussed in Sec. II C , the temporal evolution of Y 
is in addition governed by the two-particle correlations 
B and N, which, within a self-interaction like process, 
are induced through the non-pair conserving Coulomb 
couplings SH^ and 5H^ 2 '] the details of this coupling, 
Eqs. |TT|) , reflect the essential features of dielectric 
screening where: an electron (hole) polarizes its sur- 
rounding medium (described by the source terms for 
B); in turn, the hole (electron) not only feels the bare 
Coulomb term exerted by the the first carrier but also 
the field induced by the polarization cloud [cf. terms 
in the second line of Eq. in addition, this induced 

polarization cloud also affects the propagation of the car- 
rier [cf. X>-terms in the first line of Eq. (fl2|)l. Thus, the 
screening of the electron-hole interaction and the quasi- 
particle renormalizations originate from the same physi- 
cal process. 

It is interesting to compare this finding with the re- 
spective GWA and Bethe-Salpeter equation result. Let 
us first concentrate on the case of static screening, 
Eq. (|14 ) . Inspection of the quasiparticle renormalizations 
of Eq. ( 1J) with the results of the GWA approacbu reveals 
that these renormalizations closely resemble the screened 
exchange and Coulomb hole self-energy terms (note that 
this comparison is somewhat complicated because of our 
use of the electron-hole picture and the missing screening 
of the electron-hole exchange interactions). If dynamic 
screening is considered, the inelasticities w\2 of T>(u>i2) m 
Eq. ( |l2| ) are given by the differences between the light fre- 
quency u> and the electron- hole transition energies ef +eJj. 
Assuming that the exciton states are composed of elec- 
tron and hole single-particle states with energies close to 
the bandgap, one immediately notes that w « e\ + e|; 
thus, static screening of the direct Coulomb interactions 
is expected to be a good approximation. This result is 
different with respect to the GWA prescription, where the 
inelasticities of the quasiparticle renormalizations [first 
line of Eq. (|l2|)1 are the quasiparticle energies themselves, 
and the electmn-hole interaction is screened at the opti- 
cal frequency!], u) (see also Ref. for a discussion of the 
compensation between effects due to dynamical screening 
and higher-order vertex corrections; the consequences of 
such interplay within our present density-matrix descrip- 
tion will be discussed elsewhere). 

Next, we discuss the screening of the electron- hole ex- 
change interactions v ch,hc . As we saw in the discussion 
of the screening of the direct Coulomb terms, dielectric 
screening originates from a process where: the interband 
polarizations Y induce the electron-hole pair correlations 



B and N; these induced polarization fluctuations move 
independently of each other (random-phase approxima- 
tion); and finally couple back to the propagation dynam- 
ics of Y. These self-interaction-like processes result in 
the terms ©11,22) Eq. (|l3|), where the vertices 22 and 11, 
respectively, reflect the coupling of Y to the electron-hole 
pair coherence B and the back action of B, N on Y; fi- 
nally, the propagation dynamics of B and N is hidden in 
T>, which has to be determined from the Dyson-like equa- 
tion for the screened Coulomb potential W. Let us now 
return to the full expression of Eqs. (10 ll]) which, within 
the approximations adopted, contain all possible source 
terms. One immediately recognizes that neither these 
source terms nor the coupling of B, N to Y, Eq. (||), can 
reproduce a screening term of the form X> ch ' hc . Hence, 
within our approximation scheme (RPA) the electron- 
hole exchange interactions must remain unscreened. Fu- 
ture work will address possible screening contributions 
and vertex corrections beyond the random-phase approx- 
imation. 

Finally, we briefly comment on one of the well-known 
shortcomings of ou r density-matrix description. From 
Eqs. ( |l2"| ) and (B5) one readily observes that the un- 
renormalized single-particle energies e° and e h enter into 
the calculation of the screened Coulomb potential and of 
the polarization function P. This result differs from the 
corresponding GWA and Bethe-Salpeter equation result 
and is due to our neglect of higher-order correlation func- 
tions (i.e., correlations between three at, more electron- 
hole pairs). However, it can be shownfij that inclusion 
of certain types of Coulomb interactions at this level of 
many-particle correlations indeed results in a renormal- 
ization of the single-particle energies e c and e h . For the 
sake of brevity, here we shall not present the details of 
such analysis. 

When defining the Coulomb interactions of our start- 
ing Hamiltonian, Eq. (Q), we assumed that exchange in- 
teractions are fully accounted for within the DFT calcu- 
lations. Apparently, this assumption no longer holds for 
DFT calculations based on the celebrated local-density 
approximation (LDA)Jd. Inspection of Eq. (||) reveals 
that in this case the substracted exchange contributions 
would be of the form v x (r,r') = 5(r — r')v^ DA (r); as 
consequence, in SH only the Hartree contributions with 
the filled valence band would be cancelled exactly (with 
no corresponding cancellation for the exchange terms) 
and the local exchange potential would give rise to terms 
proportional to c^<$ and dc. In the equations of motion 
for the dynamic variables these local-exchange contribu- 
tions would give rise to additional couplings between vari- 
ables containing different numbers of electron-hole pairs. 
Within the spirit of perturbation theory, such couplings 
are due to the fact that the semiconductor groundstate is 
not given by simply filling the DFT-LDA states; rather, 
the non-local exchange interactions give rise to admix- 
tures of excited DFT-LDA states, which, within the 
electron-hole picture, correspond to dynamic variables 
accounting for different numbers of electron-hole pairs. 
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However, it is well known that in many cases the single- 
particle renormalizations to the DFT-LDA states simply 
lead to a rigid shift of the calculated bandgap (scissor 
operator) and that the calculated wavefunctions remain 
unchanged. In this case, the results of our analysis re- 
main valid and only the electron-hole exchange interac- 
tions V eh ' he , which correct for the missing Coulomb ex- 
change interactions because of optical excitation, should 
be computed within the local-density approximation. 

V. CONCLUSIONS AND SUMMARY 

In conclusion, starting from an ab-initio bandstruc- 
ture calculation we have shown how to describe dielec- 
tric screening within a density-matrix description. Wc 
have identified those dynamic variables and approxima- 
tions which are required at the lowest level of approxima- 
tion (random-phase approximation). We have discussed 
that static screening is expected to be a good approxima- 
tion (is contrast to the corresponding Green's function 
resultsa which seem to favour screening at the optical 
frequencies). Hence, our analysis suggests calculation of 
the optical absorption spectra feftW the solutions of the 
"excitonic eigenvalue problem" :Btffi 

(15) 

with E x and ty x the exciton energies and wavefunctions, 
respectively (with proper summation over barred in- 
dices). Here, we have assumed that single-particle renor- 
malizations (scissor operator, GW corrections) lead to a 
rigid shift of the calculated bandgap, with e^' 11 the renor- 
malized single-particle energies [see Eq. (|l2]) for the more 
complete expression]; the W's are the Coulomb interac- 
tions screened within RPA, where the dielectric function 



is computed at zero frequency (static-screening approx- 
imation) using the renormalizcd single-particle energies 
Ej 6 ' ; finally, V ch,hc are the unscreened electron- hole ex- 
change interactions; within LDA, these terms should be 
computed within the local approximation. In the optical 
spectra, the optical transition energies are provided by 
E x , and the respective oscillator strengths are given by 

IEl2*12M^| 2 . 

It is hoped that our approach might serve as an al- 
ternative to the hitherto used GWA and Bethe-Salpeter- 
equation approaches. Its advantage lies in its conceptual 
simplicity: only some basic knowledge about statistical 
mechanics and Heisenberg's equations of motion are re- 
quired (as opposed to the more sophisticated framework 
of Green's functions) . Finally, this work might help to ex- 
tend the variety of techniques developed recently within 
the field of semiconductor optics for the description of 
the nonlinear optical response to first-principles studies. 
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APPENDIX A 

In this Appendix we discuss the transformation of 
Eq. (||) to the electron-hole picture. The relation be- 
tween the field operators xp and c, d is given by: 



^(r)=^Ur)rc\+^2^(r)d 2 (Al) 

1 2 

Inserting this expression into Eq. (||), we obtain after 
some straightforward algebra: 



AffW - V oh ' cc c^Vc - V eh '- h JJdd +V oc,hc c t cdc -V hh ' hc d f ddc (A2b) 

on — V 12j g 3 o 1 «2<-gi-3 v 12,33 c l u 2 u 3 a 3 ^ v 33,21 L 3 L 3 U 2 L 1 v 33,21 U 3 U 3 U 2 L 1 \±^4U) 

6HW = i(V^ 4 h cI444 + V^ 4 h 3 C d 2 d 4C3Cl ) (A2c) 



(we assume an implicit summation over all single-particle indices); the commutation relations between c, d and the 
various Coulomb terms are: 
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[d 2 ,6HM] 
[d 2 ,SH ( 
[d 2 ,6H^] 



r ee,ee J 
11,33 3 



3 U 1 



22,33 3 3^2 



TO 

-(v. 



/t rCCCh 

V V 11,34 



V 31,14- ) 



hh _ 

33 
hh,cc 
22,33 

"\d\c 



V 



ih,h 
13,3 
he, oh 
23,32 



I) 4^3 C 1 



V 12,33 a 2 a 3 a 3 



V 



/-irhh,ch -irhh,eh^ 
_ 22.34 _ 42 



22,34 
ch,ch _f 



42,32 I 



!.'. C 1 C 3 C 3 



V 12,34 L 



3 a 4 a 2 



_ v ch,eh t J 
V 12,34 L 3 a 4 



c l 



V ee,he j 

V 11,43 "4 c 3 C l 
-trhh,hc j 7 

V 22,43 fl 4 C 3 a 2 



(A3a) 
(A3b) 
(A3c) 
(A3d) 
(A3e) 
(A3f) 



r 



Finally, the Coulomb matrix elements in Eq. (A2) are 
defined as: 



V 11.22 



d ( rr ') T^T^li^M^r'), (A4) 



sponse we consider a monofrequent laser excitation of the 
form £ e~ M *. Since Y, B, and N are driven by the light 
field, their time dependence is also e~ %u}t . Hence, keeping 
only self-interaction like processes and direct Coulomb 
terms, we find for B and N the Dyson-like equations: 



with 



d>g(r) 
*i2(r) 



[0!(r)]VI(r) 
[$(r)]*#(r) 

[#(r)^(r)]* 
^(r)^f(r). 



(A5a) 
(A5b) 
(A5c) 
(A5d) 



12,34 



-^12, 34 



V 34,43 -°12,34 



34 



Wl2 - £3 
V 43,43 -°12,34 



l + i0+ 

hc,ch 7V7- 
43,34 iV 12,34 



U>12 



e\ + i0+ 



2.31 



(Bl) 



APPENDIX B 



In this Appendix we sketch the derivation of our final 
expression (|l2|) . Without loss of generality, in linear re- 



with uj\2 = lu — e? 



e 2 and 



J 



D 12,34 



V 11,34 r 12 



-v j-hh,ch 
V 22,34 



n 2 



W12 



j0+ 



(B2) 



Here, the infinitesimal imaginary part z0 + has been introduced to ensure causality. Next, it turns out to be convenient 
to introduce the mixed representations -Bi2(t) = YIm ^34( r ) -^12,34 and Ni2(r) — J2 34 ^4l( r ) -^12,34- The coupling 
between B, N and Y [third and fourth line in Eq. (|9|)] can then be cast to the form: 



d(rr') 



\r — r 



Tl (*Vt(r)[B l3 (r') + N l2 (r')] - *g(r)[B l3 (r') + N 12 (r')}) 



(B3) 



(note our neglect of exchange interactions). The sum of B and N, which appears on the right-hand side of Eq. flB3| ), 
is obtained by summing Eqs. (Bl) and transforming to the mixed representation for B and N. We obtain: 



B 12 (r) + N 12 (r) = f d(ss') P(r, s;u 12 )—!—{B 12 (s') + N 12 (s')} + B° 12 (r), 
J \ s - s \ 

with the polarization in random-phase approximation: 



(B4) 



P(r,r';u;) = J2 



34 



*&(r)*S§(r') 



i0+ 



} + e§ + e 1 ? + i0+ 



(B5) 
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One readily recognizes that Eq. (B4) closely resembles Dyson's equation for the screened Coulomb potential W, 
which, in shorthand notation, reads W = V + WPV. Together with: 



W(r,r';w)= / ' ds —L-K(a,r';u), 
J \r — s\ 



(B6) 



wher e K is the inverse dielectric function]^ we immediately obtain (B + N) = KB . Inserting this expression into 
Eq. (B3) we find after some straightforward calculations: 



1 



w 12 ~ e 3 — e 4 
1 



i0^ 



-V 



ch.ee . iTTiin.cn/ \ 

34JI n 2 +W 2 2; 34 (^ 1 2) 



hh.ch, 



i - c c <rh 4_ in+ 34,22 12 
-'12 ~ e 3 — £ 4 ' ZU 



e? - + iQ+ 34,22 



yCh.hh y. _ w hh,0h, X 

v 34.22 12 vv 22.34 V^12y 



1 



w 12 — e 3 — £ 4 



h +i0 + V 34,11 r 12' 



with the matrix elements of W defined in analogy to Eq. (A4). Finally in case of time-reversal symmetry the 
wavefunctions can be chosen real. Hence, the polarization P is of the forma 



(B7) 



P(r, r'; u) = / du>' D (r, r'; a/) [ T 

Jo uj — uj 



1 



+ i0+ uj + uj' + i0+- 



(B8) 



with D {r, r'\ uj) = ^ 34 ^f^(r)^^(r')8(u> — e 3 — £4), and we arrive after some straightforward algebra at Eq. (p^). 
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